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We posit that accreting compact objects, including stellar mass black holes and neutron 
stars as well as supermassive black holes, may undergo extended periods of accretion dur- 
ing which the angular momentum of the disk at large scales is misaligned with that of 
the compact object. In such a scenario, Lense-Thirring precession caused by the rotating 
compact object can dramatically affect the disk. In this presentation we describe results 
from a three-dimensional relativistic magnetohydrodynamic simulation of an MRI tur- 

■ bulent disk accreting onto a tilted rapidly rotating black hole. For this case, the disk does 
I ' not achieve the commonly described Bardeen-Petterson configuration; rather, it remains 

O ' nearly planar, undergoing a slow global precession. Accretion from the disk onto the 

hole occurs predominantly through two opposing plunging streams that start from high 
r/^ , latitudes with respect to both the black-hole and disk midplanes. This is a consequence 

■ of the non-sphericity of the gravitational spacctime of the black hole. 



1. Introduction 



Recently, our group has been working on three-dimensional, fully relativistic simu- 
lations of tilted thick disks, first in the hydrodynamic limit 1 and here in full general 
relativistic MHD (GRMHD). A general relativistic treatment ensures that all im- 
portant relativistic features, such as the cusp in the potential and the Einstein and 
Lense-Thirring precessions of the orbits, as well as any higher order couplings of 
these features, are treated properly. The inclusion of magnetic fields is important 
because it is now widely believed that local "viscous" stresses are generated by 
turbulence that results from the magnetorotational instability (MRI). 2 

This work is carried out using the Cosmos++ astrophysical magnetohydrody- 
namics code. 3 For this work, the key features of Cosmos+-|- are the flexible mesh 
structure, adaptive gridding, and three-dimensional general relativistic MHD capa- 
bilities. This simulation is carried out on a spherical-polar mesh with an equivalent 
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resolution of 64 3 zones. The full set up of this problem and higher resolution simu- 
lations will be described in a forthcoming paper (Fragile et al., in preparation). 

2. Results 

Perhaps the most striking feature of the simulated disk at late times are the two 
opposing plunging streams that start from high latitudes both with respect to the 
black- hole symmetry plane and the disk midplane. Figure [T] shows the structure of 
the full disk and a zoomed in view of the plunging region with the plunging streams. 
Note that one stream remains entirely above the black-hole symmetry plane, while 
the other remains below. As material passes through this plunging stream it under- 
goes strong differential precession. The precession totals approximately 180°; thus, 
the material enters the black hole from the opposite azimuth from which it began 
its plunge. 



Fig. 1. Logarithmic density plot from the final data dump (t = 10 orbits at the pressure maximum 
of the initial torus). This figure is oriented such that the black- hole spin axis points straight up. 
The 2-axis of the triad gives the original orientation of the disk angular momentum vector. The 
main body of the disk shows very little warping or realignment toward the symmetry plane of 
the black hole. The plunging region (shown in inset) shows two opposing, high-latitude streams of 
material connecting the disk to the horizon (indicated by arrows). 

The interesting question from Figure [T] is why do these opposing accretion 
streams start from such high latitude with respect to the black-hole symmetry 
plane and disk midplane? Fundamentally this is due to the aspherical nature of the 
gravitational spacetime of the rotating black hole. A nice way to illustrate this is 
to consider the dependence of riscOi the coordinate radius of the innermost stable 
circular orbit, on inclination. Figure [2] illustrates this dependence for a few selected 
cases of a. The key point of the plot is that orbital stability around a rotating black 
hole is strongly dependent on the inclination of the orbit. Notice that the unstable 
region increases monotonically for increasing inclination. 
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Fig. 2. Plot of the inclination dependence of 
risco f° r black-hole spins a = 0, 0.5, 0.9, and 
0.998. Inclinations < i < 90° represent pro- 
grade orbits, whereas inclinations 90° < i < 
180° represent retrograde orbits. 



Fig. 3. Meridional plot (tp = 0) through the 
final dump of the simulation showing the log- 
arithm of density overlaid with a polar plot of 
the "ISCO surface" for an a = 0.9 black hole. 
Notice that the plunging streams start near 
the largest cylindrical radius (r cos $) of this 
surface (indicated by arrows). 



We can make better use of the information in Figure[5]by converting it to a polar 
plot and overlaying it on a plot of data from the simulation, as is done in Figure [3J 
The polar plot creates an "ISCO surface" , which gives a clear indication of where the 
most unstable regions are. Note that the plunging orbits highlighted previously start 
near where the disk first encounters the ISCO surface. More precisely, the streams 
start near the largest cylindrical radius (rcosi?) of the ISCO surface, measured 
with respect to the angular momentum axis of the disk. The increase in risco 
with inclination then explains why these plunging orbits start at high inclinations 
relative to the black-hole symmetry plane and the disk midplanc, since this is where 
the gravitational potential is most unstable. 
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